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A supersymmetric path integral representation is devel- 
oped for stochastic processes whose Langevin equation con- 
tains any number A'' of time derivatives, thus generalizing the 
Langevin equation with inertia studied by Kramers, where 
N — 2. The supersymmetric action contains A'^ fermion fields 
with first-order time derivatives whose path integral is evalu- 
ated for fermionless asymptotic states. 

1. For a stochastic time-dependent variable Xt obeying 
a first-order Langevin equation 



Lt[x] = it + F{xt) = T]t, 



(1) 



driven by a white noise rjt with {rjt) = 0, {rjtrjt') = ^u' , 
the correlation functions {xt-^ ■ ■ ■ xt^) can be derived from 
a generating functional 



Z[J] = {e'l''"^} 



VxAe 



-Sb~\-i J dtjx 



(2) 



with an action Sb — ^JdtL^, and a Jacobian A — 
det Sx^i Lt ■ We denote by S^^, Lt the functional deriva- 
tive of Lt[x] with respect to its argument. Explicitly: 



The supersymmetry implies QS = 0. 

The determinant (^) should not be confused with the 
partition function for fermions governed by the Hamil- 
tonian associated with the action (^. Instead of a 
trace over external states it contains only the vacuum- 
to-vacuum transition amplitude for the imaginary-time 
interval under consideration. In the coherent state rep- 
resentation, ct and Ct are set to zero at the initial and 
final times, respectively 

The path integral (H) can also be rewritten in a canon- 
ical Hamiltonian form by introducing an auxiliary Gaus- 
sian integral over momentum variables pt, and replac- 
ing Sb by — J dt{p1/2 - iptLt). The generator 
of supersymmetry for the canonical action is = 
Jdt(ict6xt +PtSct)- This form has an important advan- 
tage to be used later that it does not depend explicitly 
on Dt, so that the above analysis remains valid also for 
more general colored noises with an arbitrary correlation 
function {rjat-qbt') = {Dab)tt' 7^ 5tt' ■ 

Inserting (|^) into (H), the generating functional be- 
comes 



5x^,Lt=[dt + F'{xt)\5tt'- The time variable is written as Z[J]^ {e ^ '^^'^'') = j VpVxVcDc e~ 
a subscript to have room for functional arguments after J 



-S f+i J dtJx 



a subscript to have room for functional arguments after 
a symbol. It was pointed out by Parisi and Sourlas |^] 
that by expressing the Jacobian A as a path integral over 
Grassmann variables 



(6) 



A = / VcDce' 



(3) 



This representation makes the stochastic process (|l|) 
equivalent to by supersymmetric quantum mechanical 
system in imaginary time. In the supersymmetric formu- 
lation of a stochastic process, there exists an infinity of 
Ward identities between the correlation functions which 
can be collected in the functional relation 



with a fermionic action 



5/ = y dtdt'ct5x^,LtCt. = j dtct{dt+F')ct , (4) 

the combined action S = Sb+Sf becomes invariant under 
supersymmetry transformations generated by the nilpo- 
tent {Q^ — 0) operator 



Q = j dt (icSx — iLtSc) 



(5) 
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valid for an arbitrary functional $ = ^[p,x,c,c]. The 
Ward identities simplify a perturbative computation of 
the correlation functions. 

A proof of the equivalence of (0) to (||) requires a regu- 
larization of the path integral, most simply by time slic- 
ing. This is not unique, since there are many ways of dis- 
cretizing the Langevin equation (|l|). If one sets — ie, 
for i = 0, 1, 2, . . . , M, Xi = Xt^ , and Fi = F(xi), then the 
velocity x may be approximated by {xi — Xi^i)/€. On 
the sliced time axis, the force F(xt) may act at any time 
within the slice (ti,ti-i), which is accounted for by a pa- 
rameter a and a discretization F aFi + (1 — a)Fi^i. 
Note that the discretized Langevin equation is assumed 
to be causal, meaning that given the initial value of 
the stochastic variable xq and the noise configurations 
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rio,r]i, ...,r]M-i, the Langevin equation uniquely deter- 
mines the configurations of the stochastic variable at later 
time, xi, X2, xm ■ The simplest choice of the right- 
hand side of the Langevin equation compatible with the 
causality is to set it equal to 77^-1. In general, one can 
replace rji^i by Aj-i.i_i77i_i with A being an or- 
thogonal matrix, A^A = 1. The latter is just an evi- 
dence of the symmetry of the stochastic process with the 
white noise with respect to orthogonal transformations 

Some specific values of the interpretation parameter 
a have been favored in the literature, with a = or 
1/2 corresponding to the so-called Ito- or Stratonovich- 
related interpretation of the stochastic process (|l|), re- 
spectively In the time-sliced path integral, these 
values correspond to a prepoint or midpoint sliced ac- 
tion H,^. Emphasizing the a-dependence of the sliced 
action, we shall denote it by . This action is super- 
symmetric for any a: = 0, and the sliced genera- 
tor = ^iiicidxi +Pidci) turns out to be independent 
of both the interpretation parameter a and the width e 
of time slicing Q . A shift of a changes the action by the 
Q-exact term, 

S^+Sa = S^ + 5aQ"G, (8) 

where G is a function of a and a functional of p, x, c, c. 
This makes the Ward identities independent of a, i.e. 
on the interpretation of the Langevin equation. In- 
deed, setting 5a — —a we find e~"^» = e~^« e"*^ * = 
g-So (^i _|_ QHf^'^y Substituting this relation into (0) we 
observe that the a-dependence drops out from the Ward 
identities because of the supersymmetry Sq = and 
the nilpotency {Q")^ = 0. 

The simplest situation arises for the Ito choice, a — 0. 
Then the sliced fermion determinant A becomes a triv- 
ial constant independent of x. In the continuum limit 
of the path integral, however, this choice is inconvenient 
since then the limiting action Sq cannot be treated as an 
ordinary time integral over the continuum Lagrangian. 
Instead, Sa=o goes over into a so-called Ito stochastic in- 
tegral The Ito integral calculus ^ differs in several 
respects form the ordinary one, most prominently by the 
property J dx ^ J dtx. This difficulty is avoided taking 
the Stratonovich value a — 1/2, for which the contin- 
uum limit of 5*1/2 is ordinary integral Splitting 
as Sa — Si/2 + (a — 1/2)QG, the non-Stratonovich 
part vanishes in the continuum limit because Q does not 
depend on the slicing parameter e, whereas G is pro- 
portional to e ^ For a = 1/2, formula (H) has 
a conventional continuous interpretation as a sum over 
paths, and can be treated by standard rules of path in- 
tegration (e.g., perturbation expansion around Gaussian 
measures). The price for this is the additional fermion 
interaction, which possesses as an attractive feature the 
additional supersymmetry. 

The aim of our work is to extend this supersymmetric 
path integral representation to stochastic processes with 



higher time derivatives 

Lt = -i{dt)±t + F[xt) = m , (9) 

where 7 is a polynomial of any order — 1, thus pro- 
ducing N time derivatives on xt- This Langevin equa- 
tion may account for inertia via a term mdt in j{dt), 
and/or an arbitrary nonlocal friction J dr^rit-T — 
Yln=o ^nd^^xt where in = j dTjr{-T)"- /n\. The main 
problem is to find a proper representation of the more 
complicated determinant A = det S^^, Lt in terms of 
Grassmann variables. The standard formula (||), though 
formally applicable, does not provide a proper represen- 
tation of the determinant of an operator with higher- 
order derivatives because of the boundary condition 
problem. This problem is usually resolved via an op- 
erator representation of the associated fermionic system. 
In the stochastic context it has so far been discussed only 
for the single time derivative In the first-order for- 
malism, the fermion path integral can be defined in terms 
of coherent states with the above discussed vacuum- 
to-vacuum boundary conditions. Higher-derivative theo- 
ries, however, have many unphysical features, in partic- 
ular states with negative norms and it is a priori 
unclear how to define the boundary conditions for the as- 
sociated fermionic path integral. In gauge theories, the 
Faddeev-Popov ghosts give an example of a fermionic the- 
ory with higher-(second-)order derivatives. There, un- 
physical consequences of the negative norms of the ghost 
states are avoided by imposing the so called BRST invari- 
ant boundary conditions upon the path integral. For the 
above stochastic determinant with higher-order deriva- 
tives, the correct boundary condition are unknown. 

2. The solution proposed by us in this work is best 
illustrated by first treating Kramers' process where one 
more time derivative is present, accounting for particle 
inertia, i.e. where j{dt) — dt + j for a unit mass m = 1. 
Omitting the time subscript of the stochastic variables, 
for brevity, we replace the stochastic differential equation 
(^) by two coupled first-order equations 

Ly ^ V + -fV + F{X) ^ Vy, (10) 

L^^x-v^iy^, (11) 

There are now two independent noise variables, which 
fiuctuate according to the path integral 

{F[x,v]) = J VvxVuyF[x,v] 

A parameter a regulates the average size of deviations of 
X from V in Eq. (|ll|) . If we regard the basic noise correla- 
tion fmictions as functional matrices {Dn)u' — {i^nt^nt') 
for n = x,v, which act on functions of time as linear 
operators Dnft = J dt' {Dn)tt'ft', the noise correlation 
functions associated with (O) are 
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A, 



1 - cr , 



{-dt^e^^'d^^) e-^* . (13) 



Substituting (|l l| ) into (10) we find the two-derivative ver- 
sion of (^, a; + 7i + F{x) — rja, driven by the combined 
noise 



f]^ = l^v + i'x + IVx- 

This noise is white for any choice of a: 

iVcrt) = , {rj^t-qat') = Stt' 



(14) 



(15) 



Let xt [rj] be a solution of the original Langevin equation 
(^, and Xt a solution of the system ([ll|), (|lO|). The 
property ( |l5[) implies that xt[rj„] has the same correlation 
functions as Xt[rj\, for any a, thus describing the same 
stochastic process. The freedom in choosing a will later 
be used to make the effective supersymmetric action local 
in time. 

Once we have transformed Kramers' process into a sys- 
tem of coupled first-order Langevin equations ( [Til ) and 
( |lo| ) which is a trivial extension of the first-order equa- 
tion (^ to a matrix form, there obviously exists a path 
integral representation analogous to (^. It is for this 
reason that we have introduced a two noise variable and 
a fluctuating relation between x and v in Eq. ( ITTl ) . There 
is a complication though in that the noise Vx is no longer 
white since Dx is nonlocal in time. However, as ob- 
served above this does not affect the supersymmetry in 
the canonical form (^ of the path integral since the su- 
persymmetry generator does not depend on Dx (in 
contrast to Q). 

Thus, having established the supersymmetric path in- 
tegral representation of the equivalent first-order stochas- 
tic system, our strategy is to integrate out all auxiliary 
variables we have introduced and, thereby, derive the 
proper boundary conditions for the fermionic path in- 
tegral in the higher order stochastic processes as well as 
to construct the supersymmetry generator in the initial 
configuration space. 

To prepare the notation for the later generalization 
to a stochastic differential equation with N derivatives, 
we rename the variables x and v as Xm with a = 1,N^ 
and for the moment N = 2. Only the equation for xjv 
contains the force F = F{xi). The other equation just 
establishes a fiuctuating equality between xi and X2, the 
original process being described by x = xi. Inserting 
the stochastic equations (|l^) and (|l^) into the exponent 
of (12), we repeat the previous procedure and, choosing 
midpoint slicing with a = 1/2 a la Stratonovich, we ob- 
tain the path integral representation of the generating 
functional 



Z[J]= JvpVxVcVce'^"^'!'^"'' , 

•S-^ ^ Jdt (^^PnDnPn " iPnLn + Cn5x^Ln(- 



(16) 



n=v' 

It is readily verified that = 0, using the fact that 

Eferrm CmCk{SzJz„Lm)cn ~ 4 = due to the Grass- 
mann nature of c„. Explicitly, the Fermi part of the 
action reads 

Sf= Jdt [CxCx+CvCv - CxCv+CyCxF' {x)+^CvCv] ■ (18) 

The Gaussian path integral over momenta in ([l^) has 
a meaning without time slicing, and can be performed 
to recover the Lagrangian version of the supersymmetric 
action 
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1 



(19) 



The associated generator of supersymmetry is obtained 
from (|l^) by substituting into the solutions of the 
Hamilton equations of motion p„ = iD^^L^ which ex- 
tremize {Sp^S^ = 0), leading to 



Q 



2 

E 

n=l 



dt{iCnSx„ - iD^^LnSc^) 



(20) 



The final step consists in integrating out the auxiliary 
variable X2 = v, which only ap pears quadratically in the 
bosonic part of the action. Making use of the explicit 
form of Dn given in (p^, we obtain the Lagrangian form 
of the supersymmetric action 

= jdt ^Lt (^1 + ^^dtDxdt^ Lt + Sf , (21) 

where Lt is now the left-hand side of the initial equation 
(^), for the Kramers process at hand: Lt — Xt + jxt + 
F{xt). At this stage, the effective action is nonlocal in 
time. Now we take advantage of the freedom in choosing 
the parameter a. We go to the limit cr —> 0, in which 
case Dx ~ cr vanishes, reducing the action to the local 
form 



So - / dt^L^, + Sf 



(22) 



To find the generator of supersymmetry in this repre- 
sentation, we omit Sx^/ = Sy in (20), and replace xn = v 
by the solution of the equation of motion 



5vS 



-D-^L, 



1-a 



i-dt + i)Ly = 



(23) 



The generator of supersymmetry is 



In the limit cr ^ 0, D^^ ^ diverges, leading to an 
exact equality Lx ~ v — x = Q, rather than the fluctu- 
ating one (pi]). To take the limit cr ^ in the operator 
( p0| ) , one must flrst substitute ( p3| ) into the would-be sin- 
gular term D~^Lx in Q, and then take the limit. The 
supersymmetry generator assumes the flnal form 
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= dt [iCxSx ~ i{~dt + l)Lt5c^ - iLfSsJ 



(24) 



The action (g2|) provides us with the desired supersym- 
metric description of processes with second-order time 
derivatives. An important feature of the supersymme- 
try generated by Q is that the supermultiplet contains 
one boson field and two fermion fields. The reason for 
this is, of course, that a boson field with N time deriva- 
tives in the action carries N particles, each of which must 
have a supersymmetric fermionic partner. The fermion 
degrees of freedom have the conventional first order ac- 
tion, which permits us to impose the vacuum-to- vacuum 
boundary conditions within the coherent state represen- 
tation of fermionic path integrals . The boundary con- 
ditions for the bosonic path integral are the causal ones: 
Xt=o = xo and Xt=o = vq = xq. 

We have circumvented the problem of the boundary 
condition for the determinant of a higher-order operator 
by enlarging the number of Fermi fields, thereby reducing 
the problem to the known one for the determinant of the 
single-derivative operator. What happens if we integrate 
out the auxiliary Grassmann variables c„, c„. In these 
variables, the action (|l^) is harmonic, driven by external 
forces Cx and CxF'{x). After a quadratic completion the 
integration with the vacuum-to- vacuum boundary condi- 
tion yields det(9t +7)- The effective action for the other 
fermion pair becomes non-local 



Sf= dt [cxCx + cx{dt + -ry\F'{x)cx)] 



(25) 



The total effective action S = Sb + S f is still supersym- 
metric. The supersymmetry is generated by the operator 
(p^), if the last term in Q is dropped. The action ( |25| ) 
is the first-order action. So with the vacuum-to-vacuum 
boundary condition the integral over Cx,Cx would also 
give a determinant. Thus we get the representation 



det[dt + 7] det[at + {d + j)-^F'{x)] 



(26) 



Invoking the formula for the determinant of a block ma- 
trix, the non- locality in the second determinant can be re- 
moved, while maintaining the linearity in the time deriva- 
tive 



det 



a* + 7 F' 

-I dt 



(27) 



which is exactly the determinant arising from the two- 
noise process ( |To| ) , (jl^) . In this way we have represented 
the determinant of the second-order operator as a deter- 
minant of a first-order operator acting upon a higher- 
dimensional space for which the boundary conditions are 
known. 

Thus, with the help of two coupled equations driven 
by auxiliary noises we have succeeded in giving a unique 
meaning to the path integral representation of the 
Kramers process. The final path integral can be time- 
sliced in any desired way (prepoint, postpoint, midpoint. 



or any combination of these) — as long as the slicing is 
done equally in the bosonic and the fermionic actions. In 
Section 4, the procedure will be generalized to a friction 
coefficient 7 which is a function of x. 

3. We now generalize our construction to stochastic 
processes of an arbitrary order N. As a result we shall 
arrive at a supersymmetric extension of general higher 
order Lagrangian systems with a supermultiplet of N 
fermion fields which all possess a good quantum theory 
due to their first-order dynamics. 

Consider a system of coupled stochastic processes 



N 



Ln = XN + ^7n-ia;„ -I- F{xi) = i^n 



(28) 



n=l 



Ln ^ in - Xn+1 ^ Vn , n = N - 1 , N - 2 , . . . , 1 , (29) 

where xi = x. This stochastic process is equivalent to 
the original one if we assume the noise average as being 
taken with the weight e~^" , generalizing that in ( [l^ ) to 



S -1 



dt 



1 



N-l 



1 



V ^2 



, (30) 



where a = Y.n=i ^" ^^'^ ^» ^ YZ-n iN-md" jn = 
1. As for N = 2, equations (|28|) and ( |29| ) can b e com- 
bined into a single equation Lt = tyNt+J2n=i ^N-n^nt = 
■q^f From (|3^) follows that {ri„t) = and {ilatVa^) = Stf- 
Thus the correlation functions of the system (|2q ) are the 
same as of the original one. Note also that the combined 
noise correlation functions do not depend on the param- 
eters cr„. We shall assign some specific values to the a's 
to simplify the sequel formalism. 

The Hamiltonian path integral for the stochastic sys- 
tem ( p^ ) and (|9|) has the form ([l6|), where the label n 
runs now from 1 to N. With the same extension of the 
index sum, the operator in ( p^ ) generates supersym- 
metry. The noise correlation functions (^) are general- 



ized to Dn = ^-Y,n=i ^ri and £>„ = cr„(A 



A 



N- 



After integrating out the momenta p„ we arrive at the 
action ( p9[ ) with the extended sum, and the generator of 
supersymmetry assumes the form ( pO[ ) with the extended 
sum. 

Integrating out the auxiliary variables Xn is now tech- 
nically more involved, but the integral is still Gaussian. 
A successive integration is possible by observing that the 
fermion action does not depend on the variables Xn for 
n > 1, the stochastic process being nonlinear only in 
the physical variable xi = x. The classical equations of 
motion 5^ S = Q can be written in the form 



F>n-lLn-l 







(31) 



for n = 2, 3, 'N . Combining the equations for n = N 
and n = N — \, and the result with the equation for 
n = N — 2, and so on, we derive the relation 
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'N- 



fe=0 



(32) 



having inserted Dn — 1 — a and with n = 2, 3,..., TV. 
These expressions may be substituted into the action 
(p^), and the generator (pO|). As in the case N — 2, 
the supersymmetric Lagrangian action and the opera- 
tor Q turn out to have a smooth hmit (T„ — > 0. Since 
Dn ^ I/ctti, we see from (^2|) that L„ ^ 0, and we 
recover the physical relations i„ — Xn+i and, hence, 
9"a;. The action assumes the form (|2^), with Lt of 



Eq. 



The generator of supersymmetry becomes 



N 



dt I ciSx - 



n=l 



'N-n 

E 

. fc=0 



(33) 



For convenience, we give the fermion action explicitly: 



Sf^ dt 



N 



/ TV ^ 
\n=0 / 



(34) 



The operator ( p3| ) transforms the original stochastic vari- 
able X — xi into the Grassmann variable ci, = ici, 
whereas all the fermionic variables are transformed into 
some functions of the only bosonic variable x. The 
fermionic action (^) is constructed in such a way that 
QSf depends only on ci. The terms containing the other 
Grassmann variables are cancelled amongst each other. 
The ci-term is cancelled against the term resulting from 
QSb, i.e. Q{Sb -\- Sf) — 0. It is important to realize 
that the fermions are coupled with each other, and thus 
belong to an irreducible supermultiplet. The number of 
fermion is equal to the highest order of the time deriva- 
tive entering the bosonic action, as observed before for 
N = 2. 

4. The idea of splitting the higher order Langevin 
equation into a system of coupled first-order stochastic 
processes with a combined noise can also be applied to 
construct a supersymmetric quantum theory associated 
with the higher order stochastic process where the coeffi- 
cients 7„ are functions of Xt- We illustrate this with the 
example of Kramers' process with the friction coefficient 
being a function of the stochastic variable Xt ■ 

A straightforward replacement of 7 by 7(2;) in (|1C 
would yield a problem because the combined noise r]cr 
appears to be a function of xt, making the system (10), 
( pT| ) inequivalent to the original stochastic process (if the 
Gaussian distributions for the auxiliary noises are as- 
sumed). To resolve this problem, we take two coupled 
non-linear first-order processes 



= i} + V + \y{x) = v„ , 

= X - V + \x{x) = Vcr . 



(35) 
(36) 



The functions \x,v are subject to the condition 

a; - 7 - 1 , A„ = F - A, . (37) 

With the noise average defined b y (n2[ ) and the condition 
(^ ) , the stochastic system (|5|) , (|36|) is equivalent to the 
original system Lt = x -\- ^{x)x -\- F(x) = rj. 

The difference between ( p6| ) and (11) is just the ex- 
tra force Xx, which does not affect the derivation of the 
associated supersymmetric action. Repeating calcula- 
tions of section 2, we arrive at the supersymmetric action 
S = Sb + Sf where 



Sh 



1 



dt (x + -f{x)x + F{x)Y 



Sf= J dt [CxCx + CvCv + CxCx{^{x) - 1) -t- CyCv 

-CyCxiF'{x) - -yix) -I- 1) - CxCy] . 

The supersymmetry generator has the form 



(38) 



(39) 



Q = J dt {icx6x - iLtSc^ - i{-dt + l)Lt6cJ ■ (40) 

It is not hard to verify that QS = 0. 

If we set 7 to be independent of x in ( |39| ) , the fermionic 
action does not turn into (^8|), in contrast to what one 
might expect. The reason is that the fermionic path 
integral exhibits a large symmetry associated with gen- 
eral canonical transformations on the Grassmann phase 
space spanned by c and c. Recall that under canoni- 
cal transformations the canonical one-form Cndcn is 
invariant up to a total differential dF{c, c). Also the mea- 
sure Y[n dcndcn remains unchanged. Thus there exists in- 
finitely many equivalent supersymmetric representations 
of the same stochastic process. The situation is similar 
to the BRST symmetry |^] in gauge theories where the 
BRST charge is defined up to a general canonical trans- 
formation. This freedom can be used to simplify the 
fermionic action or the Fermi-part of the supersymmetry 
generator. 

This formal invariance of the continuum phase-space 
path integral measure with respect to canonical trans- 
formations has been studied thoroughly ||] for bosonic 
phase spaces. A regularization of the continuum phase- 
space path integral measure with respect to canonical 
transformations on a phase space which is a Grassmann 
manifold is still an open problem. 

Acknowledgment : 

The authors are grateful to Drs. Glenn Barnich and Axel 
Pelster for many useful discussions, and to Prof. John 
Klauder for comments. 



5 



[1] G. Parisi and N. Sourlas, Phys.Rev.Lett. 43, 744 (1979); 
Nucl.Phys. B206, 321 (1982); 

M.V. Feigel'man and A.M. Tsvelik, Sov.Phys. JETP, 56, 
823 (1982); Phys.Lett. 95A, 469 (1983); 
For a comprehensive review see J. Zinn- Justin, Quan- 
tum Field Theory and Critical Phenomena (2nd Edition, 
Clarendon Press, Oxford, 1993). 
[2] H. Ezawa and J. R. Klauder, Prog.Thor.Phys. 74, 104 
(1985); 

L.P. Singh and F. Steiner, Phys.Lett. 166B, 155 (1986); 
H. Nakazato, K. Okano, L. Schiilke and Y. Yamahaka, 
Nucl.Phys. B346, 611 (1990). 

[3] C.W. Gardiner, Handbook of Stochastic Methods (Springer 
Series in Synergetics, Vol. 13, Springer, Berlin, 1983). 

[4] H. Kleinert, Path Integrals in Quantum Mechanics, Statis- 
tics and Polymer Physics, World Scientific, Second Edi- 
tion, 1995. 

[5] In Section 10.5 of Ref. it is shown that the correct time 
slicing of an interaction J dtqF{q) in a path integral is of 
the midpoint type, corresponding to a = 1/2. Sometimes 
this is referred to as the midpoint prescription for defining 
the sliced action, but it can actually be derived from the 
short-time action along a classical orbit. 

[6] A. Pais and G.E. Uhlenbeck, Phys. Rev. 79, 145 (1950); 
M.V. Ostrogradsky, Mem. Acad. Sci. St-Petersburg, 6, 
385 (1850); See also; E.T. Whittaker, A Treatise on the 
Analytical Dynamics of Particles and Rigid Bodies (Cam- 
bridge University Press, Cambridge, 1959); 
and Section 17.3 in Vol. II of H. Kleinert, Gauge Fields in 
Condensed Matter, (World Scientific, Singapore, 1989). 

[7] M. Henneaux and C. Teitelboim, Quantization of Gauge 
Systems (Princeton University Press, Princeton, 1992). 

[8] J.R. Klauder, Ann. Phys. (NY), 188, 120 (1988). 



6 



